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Abstract 

We show that dielectric or even perfectly conducting cylinders can be cloaked by a uniform or 
a layered dielectric cladding, without the need of any exotic or magnetic material parameters. In 
particular, we start by presenting a simple analytical concept that can accurately describe the 
cloaking effect obtained with conical silver plates in the visible spectrum. The modeled structure 
has been originally presented in [S. A. Tretyakov, P. Alitalo, O. Luukkonen, C. R. Simovski, Phys. 
Rev. Lett., vol. 103, p. 103905, 2009], where its operation as a cloak in the optical frequencies 
was studied only numerically. We model rigorously this configuration as a multi-layer dielectric 
cover surrounding the cloaked object, with excellent agreement to the simulation results of the 
actual device. The concept of using uniform or multilayer dielectric covers, with relative permit- 
tivities larger than unity, is then successfully extended to cloaking of impenetrable objects such as 
conducting cylinders. 
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I. INTRODUCTION 



An electromagnetic cloak is a device that minimizes or even nullifies the effects of scatter- 
ing from various objects by rendering them electromagnetically invisible to a detector. One 
well-known cloaking technique is the so-called "transformation-optics" method^- which uses 
strongly inhomogeneous and highly anisotropic materials to perform complex transforma- 
tions for the incident wave. Another cloaking technique is based on the cancellation of the 
scattering effects of a magnetodielectric object with a suitable metamaterial or plasmonic 
cladding^, which neutralizes the polarization current of the primary scatterer. 

Unfortunately, the actual construction of a cloaking device is very difficult with both the 
aforementioned approaches. Suggestions to overcome such a drawback includes, e.g., the use 
of transmission-line networks or other waveguiding structures instead of fictitious materials 
with exotic properties.- With transmission-line networks, the object to be electromagneti- 
cally hidden is strongly limited in size and geometry,-^ but cloaking of impenetrable (e.g., 
perfectly metallic) objects has been shown to be possible with a set of conical conducting 
plates placed around a cylindrical cloaked region.— This metal-plate cloaking has been prac- 
tically realized in microwaves and numerically proven to be functional for a wide band of 
optical frequencies.— However, the underlying theory and the physical principles that govern 
the related phenomena have not been yet thoroughly understood. 

In this work, we introduce a new, simple, nonmagnetic cloak. The idea originates from 
an analytical model of the previously reported metal-plate cloak operating in the visible 
range.— The used model is a simple mathematical concept comprised of multiple concentric 
cylinders of different dielectric permittivities. The choice of the permittivities has been 
made by considering the tapered conical lines, under the related plane-wave excitation, as a 
series connection of capacitors. The analytically evaluated response of the described model 
is shown to have a remarkable coincidence with the corresponding results obtained from 
full-wave simulations of the actual device for various values of input parameters. 

Since the aforementioned simple model appears to resemble so successfully the real- world 
metal-plate cloaking structure, we tried to adopt an even simpler concept to cloak an impen- 
etrable, perfectly electrically conducting (PEC) volume. This is a clearly more challenging 
case since one can hide practically anything in a PEC volume by avoiding to interact with 
the environment, which does not happen for penetrable objects. We introduce uniform or 



2 



layered dielectric claddings which are optimized by sweeping the relative permittivity values 
of these covers. The scattering reduction achieved by the proposed devices is quite high, 
given their simple structure. Moreover, the demonstrated cloaking bandwidths can be con- 
sidered as wide, a property attributed to the fact that the materials required in the claddings 
are simple dielectrics with relative permittivities larger than unity. It should be noted that 
most other cloaking methods for PEC objects, discussed in the literature, require magnetic 
material properties. 



II. GENERIC MODEL OF CONCENTRIC CYLINDERS 

The proposed model to mimic the physical mechanism of wave propagation in metal-plate 
devices is fairly simple and its configuration is shown in Fig. [TJ where the used cylindrical 
coordinate system (p,<p, z) is also defined. There are just U infinitely long cylindrical layers 
each of which is assigned to an integer number u = 1, • • • , U, occupying the region r u _i < 
p < r u and being filled with magnetically inert dielectric materials of relative permittivities 
e r>u . One can clearly note that ro = b and ru = a, while the background medium is vacuum 
(with eo, po intrinsic electromagnetic parameters) and is taken as the region (p > b) of our 
configuration with e rj o = 1. The internal cylinder, which corresponds to the region (U + 1) 
(p < a), can be either penetrable with relative dielectric constant e r m+x) or PEC. The 
assumed harmonic time dependence is of the form e +J ' 27r ** (/ is the operating frequency) and 
is suppressed throughout the analysis. 

We consider a plane wave excitation of electric fielcU^: 

+ 0O 

E , mc = z e ->*o'«»* = z r n Uk p)e^, (1) 

n=— oo 

where ko = 27r/' v /eo7^o is the free-space wavenumber and J n is the n-th ordered Bessel 
function. The total field in each layer u — 0, • • • , (U + 1) possesses only a z component E u 
whose expression is given by: 

+oo 

E u = [A n , u Jn(k u p) + B ntU H n (k u p)}e in t (2) 

n=— oo 

where k u = fcoy^, H n is the n-th ordered Hankel function of the second type and A UjU , B n ^ u 
are sequences of unknown complex coefficients. After imposing the necessary boundary 
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FIG. 1: The geometry of the generic model comprised of several concentric cylindrical layers around 
the cloaked object. The cloaked object is situated in region (U + 1). 
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(4) 



conditions at p — r u , u — 0, • • • , (U — 1), one obtains the following relations connecting the 
coefficients of two adjacent layers: 



(3) 



for integer n. The explicit form of the transfer matrix—^ T n>n is shown in the two-column 
equation fll]). 

By successive application of (EJ) for u = 0, • • • , (U — 1), one arrives to: 
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Depending on what is the filling material of the core region (U + 1), namely dielectric (diel.) 
with permittivity e r rfy + i) or PEC, the following expression for the coefficients of the U region 
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(boundary conditions at p = rjj 



a) is formulated: 
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H n {k v ru) 



A n ,u ,PEC 



By inspection of ([I]), it is directly obtained that A n0 = j~ n , while the physical demand 
for bounded field into the cloaked region, is translated into: B n ^u+i) = 0. Therefore, 
the scattering field by the device into the vacuum background area can be readily derived 
through the related coefficient 5„ i0 by combining (J5]), (EJ) in each case: 



FIG. 2: (Color online) The physical configuration of a single conical metal-plate cell. One half (cut 
along the xz-plane) of a single cell is shown. 

III. OPTICAL CLOAK MADE OF CONICAL PLATES 

In this section, we are going to test the model analyzed above on how accurately does it 
describe the operation of an actual device. It has been shown that the so-called metal-plate 
cloak configuration, already constructed for radio frequencies, can work in the optical band 
as well.— The optical device is comprised of periodically stacked, solid, conical, silver plates 
(of outer diameter 2b) positioned around the cylindrical cloaked region (of diameter 2a) 
with a small air gap (of thickness g) in between. One period of this structure is shown in 





The notation [D] is used for the (v,w) element of the matrix D. 
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Fig. [2] and it is a (hollow) waveguide, with linearly decreasing height from H (at p = b) to 
h (at p — a), leading the z-polarized fields around the cloaked object. In a specific case^i, 
the constructing material of the plates is silver, the (optimal) operating frequency fo = 590 
THz, while the physical dimensions are given by: b = 113 nm, a = 50 nm, g = 15 nm, 
H = 13 nm and h = 2.5 nm. It should be stressed that the cloaked region is filled with 
silver too, while the frequency-dependent permittivity of this substance e r ^u ver = e r ^u ver (f) 
is well known.— Note that the permittivity of silver is negative close to / = /o, namely: 

^[^r, silver (fo)] — —10. 




p=b p=r u P^ a+ g P =a 



FIG. 3: The side view of a single conical metal-plate cell. 

With reference to the model analyzed above, we make use of the formulas corresponding 
to a dielectric core since the region (U + 1) is solid silver; therefore, e r ^u + i) = e r ^ si i veT . In 
addition, the U-tYi layer is the air gap; thus, e r jj = 1 and ru-i — a + 9- It is apparent 
that our model does not take into account the periodic dimension variation of the actual 
configuration with respect to z, since the stack of conical plates is replaced by a structure 
of homogeneous concentric cylinders. However, the radial inhomogeneity (contrary to the 
axial one) of the cloaking device is possible to be imitated by choosing properly the dielectric 
permittivities of the cylindrical layers. But how can one compute e r>u for u = 1, • • • , (U — 1) 
in a correct and reliable way? 

In Fig. [31 one observes the side view of half the biconical cell as the shape possesses 
cylindrical symmetry. Assume that the cylindrical layer in the vicinity of the representative 
surface p = r u (containing both the corresponding vacuum aperture and the related solid 
silver plates), would be replaced by a (locally) homogeneous dielectric, with the same ex- 
ternal dimensions, of relative permittivity e r u , whose value should be estimated. From the 
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similarity of triangles, the height ( u of the representative vacuum aperture is found equal 
to: 

< . = Jr !iZ£Z£ + k »Z*_. (8) 
o — a — g o — a — g 

If the inclination of the tapered plates is relatively low (which is the case for the specific 
choice of b, a, g, H, h), then the z-polarized electric field will be almost normal to the sloping 
boundaries separating the silver from the vacuum. As a result, the arbitrary cross section 
at p = r u can be considered as a series connection of three capacitors: two with the height 
(H — C«)/2 filled with silver and one (placed in between) with height ( u filled with vacuum. 
Accordingly, this serial cluster can be replaced by a single capacitor with a dielectric material 
of relative permittivity e r , u , based on the well-known equivalent capacitance formula: = 
(h-( u )/2 _|_ _|_ {H-(u)/2 ^ j n ^yy S sense, the rigorous expression for e ru is hyperbolic with 

^r, silver J- €r, silver 

respect to both r u , ( u and is given as follows: 

H ^r, silver /q\ 

6r,u = H -I- (c _1V ' ^ ' 

n -|- y^r, silver *-)>>u 

We use layers of the same thickness (due to the constant inclination of the silver conical 
plates), namely: r u = b — jjz?[{b — a — g) ior u — 1, ■ ■ ■ , (U — 1). After having defined all the 
parameters of the theoretical structure (e r>u , r u ) based on the actual quantities of the real 
configuration (b, a, g, H, h, e r)S u ver ), we can develop the corresponding model and quantify 
its efficiency. It is remarkable that (for the given ranges of the input parameters) all the 
relative permittivities e r ^ u are positive and greater than unity. 

A crucial quantity for the operation and the performance of a cloaking device is the total 
scattering width of the whole cylindrical structure normalized by the corresponding width 
of the uncloaked one. The smaller is the magnitude, the better the device serves its purpose. 
The quantity is defined by: 

2 

d(f) 

— (10) 



f 

a, 



rZTT B -n jnq 

JO l^n=-N D n,QJ V 



where iV is a sufficiently large integer to achieve convergence for the series and B' n denote 
the coefficients of the scattering field for the uncloaked cylinder. In Fig. H] we show the varia- 
tion of (Jnorm with respect to the operating frequency for two alternative heights h = 2.5 and 
h = 5 nm. In each case, we perform a numerical simulation of the actual metal-plate device 



with ANSYS HFSS full-wave software.— The obtained numerical results are compared to 
those derived through implementation of the corresponding multi-layered analytical model, 
in the way described above. There is an excellent agreement between the two sets of data 
despite the fact that they describe two completely different structures (simulation of an 
axially inhomogeneous real device and rigorous solution to an axially homogeneous mul- 
tilayered configuration). This remarkable coincidence between so dissimilar configurations 
demonstrates the success of the adopted model. When it comes to the results themselves, 
there are large frequency bands where the magnitude o norm takes values much smaller than 
unity. In addition, the metal-plate cloaking for h = 5 nm is functional at larger frequencies 
than the device with h = 2.5 nm does. 




operational frequency /(THz) 



FIG. 4: (Color online) The total scattering width of the metal-plate structure normalized by the 
corresponding quantity of the uncloaked object as function of the operating frequency. The HFSS 
simulations for the actual device are compared with the results obtained through the analytical 
model implementation. Plot parameters: b = 113 nm, a = 50 nm, g = 15 nm, H = 13 nm, U = 8. 

In the considered numerical example, the electrical dimension of the structure is relatively 
small which could make one think that just the omnidirectional (N = 0) term of the sums in 
f fTUj) . is sufficient to evaluate the quantity (J norm- If such an argument was correct, then the 
cloaking behavior demonstrated above would be attributed to the well-known "scattering 
cancellation" phenomenon.- 1 ^ We believe that this is not the case since the electrical size 
of the structure is not so small. To support this statement, we assume the corresponding 
"scattering cancellation" device for the silver cloaked object of radius a. This is a homo- 
geneous cladding of relative permittivity e r and radius b without an air gap. In Fig. we 
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evaluate the normalized total scattering width a norm from (fit)]) as function of e r for several 
truncation limits N. We are mainly interested in the interval < e r < 5 where smaller 
&norm are observed and it is clear that for N > 2, the result does not change significantly. 
However, the first harmonic (N = 1) is necessary to be included in the computation since 
the zeroth term is minimized (and more specifically nullified) at larger e r . The nullification 
of the omnidirectional term is the "scattering cancellation" solution and it is different from 
the observed one. 
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FIG. 5: (Color online) The normalized total scattering width of a silver rod surrounded by a 
homogeneous cladding as function of the relative permittivity of the cladding. Plot parameters: 
b = 113 nm, a = 50 nm, f = 590 THz, U = 1. 



IV. CLOAKING OF CONDUCTING CYLINDERS 

As was shown in the previous section, a multilayer as well as a single layer dielectric 
cover can be used to reduce the scattering from a cylinder composed of an e r -negative ma- 
terial. The same can be of course achieved by using the well-known "scattering cancellation 
method" . However, this technique has not been used so far to hide a conducting object by 
a layered or a uniform cover made of dielectrics with e r > 1. Here we show, by using the 
analytical model described in Section II, that such a case is possible at least for conducting 
cylinders of moderate electrical size. 

The electromagnetic fields weakly penetrate the cloaked silver cylinder studied in the 
previous section, so it is expected that also conducting cylinders, i.e., impenetrable objects, 
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FIG. 6: Sweep of e r to find the optimal value for cloaking at the frequency /o with fixed values of 
a and b. The uniform (U = 1) cloak has a thickness (b — a) (with b = 2a = Ao/5) and is made of a 
dielectric material with e r . 

could also be partially cloaked with simple multilayer or even uniform covers. To present a 
generalized case, in the following we normalize the dimensions and the results to the free- 
space wavelength A . We study an electrically thin, PEC cylinder with radius a = A /10 
(the electrical size of this cylinder is of the same order as the cloaked object in the previous 
section). We assign a dielectric cover around this cylinder and consider three cases: (i) 
a uniform cover with constant e r , (ii) a multilayer cover with linearly varying e rjU , (iii) a 
multilayer cover with hyperbolically varying e T)U (as in 0). For simplicity, we choose b = 2a 
in all the studied cases. However, the choice of b and a can be made freely; the resulting 
value of (optimal) e r ^j of the cover material simply changes if the values of a and b are 
changed. 

In the first case, we assign a constant e r to the material surrounding the PEC cylinder. 
The analytical model of Section II is used to optimize the value of e r (see Fig. [6]). It is 
evident that for these values of a and b, the optimal value is e r = 5.42 and with that, the 
normalized total scattering width of the cloaked PEC cylinder is slightly less than 0.4, i.e., 
the cloak reduces the total scattering width a norm of the PEC cylinder by more than 60 %. 

The frequency dependence of the normalized total scattering width is illustrated in Fig. [7J 
demonstrating a reasonably broadband cloaking effect: the relative bandwidth where the 
total scattering width of the PEC cylinder is reduced to the half or less, is about 21 %. It is 
quite interesting to note that moderate losses do not deteriorate the cloaking effect; with a 



10 




0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 
normalized operational frequency/// Q 

FIG. 7: (Color online) Normalized total scattering width as a function of the normalized frequency. 
The uniform (U = 1) cloak has a thickness (b — a) (with b = 2a = Ao/5) and is made of a dielectric 
material with e r = 5.42. 

TABLE I: Values of relative permittivities with linearly changing e r . 



u 


1 


2 


3 


4 


5 




3.22 


5.44 


7.66 


9.88 


12.1 



loss tangent of 0.01, the cloaking effect is actually slightly improved compared to the lossless 
case. 

In the second scenario, we model the case of Fig. [1] for U = 5 and the value of e r>u 
depending on the layer number u linearly. To find the optimal values of e r ^ u for cloaking at 
/ — /o> we pl°t Pnorm as a function of e rj 5, i.e., the maximum value of e r , while e r>u (for u = 
1, 2, 3, 4) is linearly decreasing from e r $ to e r ^ = 1. The result shown in Fig. E] demonstrates 
that for the chosen dimensions a and b, the lowest positive value of e r ^, corresponding to a 
minimum in the normalized total scattering width, is e r ^ = 12.1. In that case, the relative 
permittivities of the five-layer cloak are as shown in Table I. 

With the values of Table I, the normalized total scattering width as function of the fre- 
quency looks as shown in Fig. [9j In the same graph, the curve for the constant-permittivity 
cladding is depicted for comparison and also the corresponding numerical results obtained 
with ANSYS HFSS (showing good agreement with our analytical findings) have been at- 
tached. 

Finally, we study the third scenario, i.e., the case with e r depending hyperbolically on 
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FIG. 8: Sweep of e rj u=5 to find the optimal value for cloaking at the frequency /o with fixed values 
of a and b = 2a = Aq/5. Linear and hyperbolic profiles (U = 5) are considered. 




FIG. 9: (Color online) Normalized total scattering width as a function of the normalized frequency 
for constant, linear, and hyperbolic profile (U = 5). The circles, squares, and diamonds denote the 
numerical results of each of the previous cases, respectively. 

the layer number u = 0, • • • , U, where again U = 5. The variation of the normalized total 
scattering width a norm with respect to e rj5 , is shown in Fig. [HI resulting in the value e r>5 = 128 
for optimal operation at Jq. With this value, the relative permittivities of the five-layer cloak 
are as shown in Table II. The frequency dependence of the normalized total scattering width 
is illustrated in Fig. [9j Again, we verify our analytical evaluations by plotting the simulation 
results in the same figure. 

Comparing the three curves of Fig. M, we can conclude that changing the profile of e r , u 
from constant to linear and hyperbolic, improves the obtained scattering reduction at the 
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TABLE II: Values of relative permittivities with hyperbolically changing e r . 



u 


1 


2 


3 


4 


5 




1.25 


1.66 


2.47 


4.85 


128 



frequency /o, but at the same time, the bandwidth of efficient cloaking decreases. Concerning 
practical issues, the constant and linear profiles are easily realizable (the required values of 
€ rtU are feasible), whereas the hyperbolic profile is far from practical. However, the cloaking 
performance with the constant profile is not much different from the linear profile, so it may 
not be worth the increased complexity to use even the linear one. 

It is clear that cloaking of PEC objects with simple dielectric covers is far from ideal 
cloaking such that is in theory possible with, e.g., "transformation-optics". However, it is 
evident that the cloaking efficiencies presented in this work, are comparable to the exper- 
imental and numerical results obtained with various other cloaking techniques that have 
been realized with composite metamaterials.— ^ 



V. CONCLUSIONS 



We have presented a very simple analytical concept, based on transfer matrices at mul- 
tilayered cylindrical structures. It has been found that this model describes accurately the 
previously reported cloaking effect obtained with conical silver plates in the visible frequency 
range. The effectiveness of the analytical model is verified by comparing the results of the 
normalized total scattering widths originating from it, to results obtained by numerical sim- 
ulation software. The fidelity of the proposed concept allows it to be used in device design 
and to save computational resources due to its simplicity. The same analytical model is 
also used to demonstrate that, surprisingly, cloaking of impenetrable (perfectly conducting) 
objects is possible with simple dielectric covers whose relative permittivity surpasses unity. 
Such a property renders this type of electromagnetic configurations easily realizable, unlike 
most other cloaking devices reported in the literature. 
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